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The full phonon dispersion relations of lead titanate and lead zirconate in the cubic perovskite 
structure are computed using first-principles variational density-functional perturbation theory, with 
ab initio pseudopotentials and a plane- wave basis set. Comparison with the results previously 
obtained for barium titanate shows that the change of a single constituent (Ba to Pb, Ti to Zr) has 
profound effects on the character and dispersion of unstable modes, with significant implications for 
the nature of the phase transitions and the dielectric and piezoelectric responses of the compounds. 
Examination of the interatomic force constants in real space, obtained by a transformation which 
correctly treats the long-range dipolar contribution, shows that most are strikingly similar, while 
it is the differences in a few key interactions which produce the observed changes in the phonon 
dispersions. These trends suggest the possibility of the transferability of force constants to predict 
the lattice dynamics of perovskite solid solutions. 



I. INTRODUCTION 

In the family of perovskite ABO3 compounds, a wide 
variety of distorted variants of the high-temperature cu- 
bic structure are observed as a function of composition 
and temperature jljj. First-principles density- functional 
energy-minimization methods have proved to be gener- 
ally quite accurate in the theoretical prediction of the 
ground state structure type and structural parameters of 
perovskite oxides |2| . According to the soft-mode theory 
of structural phase transitions, the ferroelectric phases 
can be related to the high temperature symmetric struc- 
ture by the freezing-in of unstable zone-center phonons. 
However, to predict finite temperature behavior at phase 
transitions, as well as the temperature dependence of the 
dielectric and piezoelectric responses of the compounds, 
it is necessary to have information about the energies of 
non-uniform instabilities and low energy distortions of 
the cubic perovskite structure ||. For this purpose, the 
calculation of the phonon dispersion relations through 
density-functional perturbation theory (DFPT) |^|,|| has 
been found to be extremely useful, allowing easy identi- 
fication of the unstable phonon branches and their dis- 
persion throughout the Brillouin zone. This information 
also permits the investigation of the geometry of localized 
instabilities in real space, which can be directly related 
to the anisotropy of the instability region in reciprocal 
space To date, phonon calculations away from 

the zone center have been reported for only a few indi- 
vidual compounds, with full phonon dispersion relations 
given for KNb0 3 @, SrTi0 3 @ and BaTi0 3 @, and 
selected eigenmodes for PbTi0 3 [M . 

In this paper, we present the first comparative study 
of the full phonon dispersion relations of three different 
cubic perovskites: BaTi0 3 , PbTi0 3 and PbZr0 3 , com- 
puted from first principles using the variational DFPT 
method as described in Section II. These compounds have 
been chosen both because of their scientific and techno- 



logical importance, and because they allow us to inves- 
tigate the effects on the lattice dynamics of substituting 
one cation with the other atoms unchanged. As we will 
see in Sections III and IV, these substitutions lead to 
very pronounced differences in the eigenvectors and dis- 
persions of the unstable phonons. The origin of these 
differences will be clarified in Section V, where, using 
a systematic method previously applied to BaTi0 3 ||, 
the reciprocal space force constant matrices are trans- 
formed to obtain real-space interatomic force constants 
for all three compounds. There, we will see that the 
differences can be attributed to trends in cation-oxygen 
interactions, with other IFCs remarkably similar among 
the three compounds. These trends and similarities, and 
their implications for the study of the lattice dynamics 
of solid solutions, are discussed in Section V. Section VI 
concludes the paper. 



II. METHOD 

The first-principles calculations for BaTi0 3 follow the 
method previously reported in Ref. M. For PbTi0 3 and 
PbZr0 3 , calculations were performed within the Kohn- 
Sham formulation of density functional theory ffu, us- 
ing the conjugate-gradients method |12| . The exchange- 
correlation energy functional was evaluated within the 
local density approximation (LDA), using the Perdew- 
Zunger parametrization [O] of the Ceperley-Alder H] 
homogeneous electron gas data j^]. The "all-electron" 
potentials were replaced by the same ab initio pseudo- 
potentials as in Refs. Jl(j] and |L6]]. The electronic wave- 
functions were expanded in plane waves up to a kinetic 
energy cutoff of 850 eV. Integrals over the Brillouin zone 
were approximated by sums on a 4 x 4 x 4 mesh of special 
fc-points [ yjj . 

The optical dielectric constant, the Born effective 
charges and the force constant matrix at selected g-points 
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of the Brillouin zone were computed within a variational 
formulation of density functional perturbation the- 
ory 0| . The phonon dispersion curves were interpolated 
following the scheme described in Ref. |l^.[nj. In this 
approach, the long-range character of the dipole-dipole 
contribution is correctly handled by first subtracting it 
from the force constant matrix in reciprocal space and 
treating it separately. The short-range contribution to 
the interatomic force constants in real space is then ob- 
tained from the remainder of the force constant matrix in 
g-space using a discrete Fourier transformation. In this 
work, the short-range contribution was computed from 
the force constant matrices on a 2 x 2 x 2 centered cu- 
bic mesh of g-points comprised of T, X, M, R and the A 
point halfway from T to R [^(J. From the resulting set 
of interatomic force constants in real space, the phonon 
spectrum can be readily obtained at any point in the 
Brillouin zone. 

Our calculations have been performed in the cubic per- 
ovskite structure. For PbTiOs, the optimized LDA lat- 
tice parameter (3.883 A) slightly underestimates the ex- 
perimental estimated value of 3.969 A |2l[] and we have 
decided to work, as for BaTiC>3 ||, at the experimental 
volume. For PbZrC>3, we have chosen, as in Ref. p2[ , 
to work at the optimized lattice parameter of 4.12 A, 
which is nearly indistinguishable from the extrapolated 
experimental value of 4.13 A |23|. 



III. DIELECTRIC PROPERTIES 

Knowledge of the Born effective charges (Z*) and the 
optical dielectric tensor (eoo) is essential for describing 
the long-range dipolar contribution to the lattice dynam- 
ics of a polar insulator. In Table [|, we present results 
for PbTiC>3 and PbZrC>3, computed using the method 
described in Section II, and for BaTiC>3, obtained pre- 
viously in Ref. ||. The effective charges have been cor- 
rected following the scheme proposed in Ref. J|4|] in order 
to satisfy the charge neutrality sum rule. Our results dif- 
fer by at most 0.09 electrons with values reported for 
cubic PbTiOs and PbZr03 usin g sligh tly different meth- 
ods and/or lattice constants mO 16.29]. 



As usual in the class of perovskite ABO3 compounds, 
the amplitudes of some elements of the effective charge 
tensors deviate substantially from the nominal value ex- 
pected in a purely ionic picture (for a review, see for in- 
stance Ref. ]|6) ) . This effect is especially pronounced for 
the Ti and associated 0|| charges in BaTiOs and PbTiOs. 
It reflects the sensitivity to atomic displacement of the 
partially-covalent character of the Ti-0 bond ]27[] . In 
contrast, the effective charge of the Zr atom in PbZr03 
and the associated Oy charge are significantly closer to 
their nominal ionic values of +4 and —2, respectively. 
The Zr effective charge is comparable to that reported 
recently for Zr02 (+5.75) p8j. In addition, in both 



PbTiOs and PbZr03, the anomalous contribution to the 
Pb charge, beyond the nominal charge of +2, is more 
than twice as large as that for the Ba charge in BaTi03. 
This feature, too, with a concomitent increase in the 
magnitude of Zq ± , reflects the sensitivity to atomic dis- 
placement of the partially covalent character of the bond 
between lead and oxygen. 

Within the LDA, the computed optical dielectric con- 
stant (Table |l|) usually overestimates the experimental 
value. The error is of the order of 20% in BaTi03 [30], 
for which the extrapolated experimental value is 5.40 [311, 
consistent with analogous comparisons in KNb03 |32f l 
and SrTi03 |J. For PbZr03, it appears that experi- 
mental data is available only for the orthorhombic phase, 
where the value is about 4.8 f33(| , significantly less than 
the value of 6.97 which we have computed for the cubic 
phase. For PbTiOs, our value of 8.24 is comparable to a 
recent first-principles result of 8.28 obtained using a dif- 
ferent method |3^| . In contrast to the other perovskites, 
this represents a slight underestimate of the extrapolated 
experimental value of 8.64 reported in Ref. pE |. 

The origin of the LDA error in the optical dielectric 
constant is a complex question. It arises at least partly 
from the lack of polarization dependence of the approx- 
imate exchange-correlation functional pjJj. In the cu- 
bic phase of the perovskite ferroelectrics, the comparison 
with experiment is also complicated by the fact that the 
high-temperature cubic phases for which the measure- 
ments are made do not have a perfect cubic perovskite 
structure, as assumed in the calculations. In fact, the ob- 
served cubic structure of BaTi03 and PbTi03 represents 
the average of large local distortions. The character of 
these distortions depends strongly on the material, which 
could well have different effects on the observed optical 
dielectric constant. 

For the BaTiOs phonon dispersion, it has been checked 
that the inaccuracy in only significantly affects the 
position of the highest longitudinal optic branch, while 
other frequencies are relatively insensitive to the ampli- 
tude of the dielectric constant Q . The effects of possible 
discrepancies for the other two compounds are likely to 
be similarly minor. 



IV. PHONON DISPERSION CURVES 

In this section, we describe phonon dispersion relations 
for BaTi0 3 , PbTi0 3 and PbZr0 3 , providing a global 
view of the quadratic-order energy surface around the 
cubic perovskite structure. The calculated phonon dis- 
persion curves along the high symmetry lines of the sim- 
ple cubic Brillouin zone are shown in Fig. ^. The un- 
stable modes, which determine the nature of the phase 
transitions and the dielectric and piezoelectric responses 
of the compounds, have imaginary frequencies. Their 
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dispersion is shown below the zero- frequency line. The 
character of these modes also has significant implications 
for the properties of the system. This character has been 
depicted in Fig. |l| by assigning a color to each eigenvalue, 
determined by the percentage of each atomic character in 
the normalized eigenvector of the dynamical matrix (red 
for A atom, green for B atom and blue for O atoms) p5[ . 

Barium titanate (BaTiOs) and potassium niobate 
(KNbOs) both undergo a transition sequence with de- 
creasing temperature through ferroelectric tetragonal, 
orthorhombic and rhombohedral (ground state) struc- 
tures, all related to the cubic perovskite structure by the 
freezing-in of a polar mode at T. The main features of 
the phonon dispersion of BaTi03 have been previously 
discussed in Ref. § and are very similar to those of 
KNbC>3 ||. The most unstable mode is at T, and this 
mode, dominated by the Ti displacement against the oxy- 
gens (Table H), is the one that freezes in to give the ferro- 
electric phases. However, the instability is not restricted 
to the r point. Branches of Ti-dominated unstable modes 
extend over much of the Brillouin zone. The flat disper- 
sions of the unstable transverse optic mode towards X 
and M, combined with its rapid stiffening towards R, con- 
fine the instability to three quasi-two-dimensional "slabs" 
of reciprocal space intersecting at T. This is the finger- 
print of a "chain-like" unstable localized distortion for 
the Ti displacements in real space . Except for these 
modes, all the other phonons are stable in BaTiC>3, which 
makes the behavior of the unstable branches relatively 
easy to understand. 

Lead titanate (PbTiOs) has a single transition to a 
low-temperature ferroelectric tetragonal structure, re- 
lated to the cubic perovskite structure by the freezing-in 
of a polar mode at T. The phonon dispersion of PbTiC>3 
shows similar features to that of BaTiC>3, with some im- 
portant differences. As in BaTiC>3, the most unstable 
mode is at T, consistent with the observed ground state 
structure. However, the eigenvector is no longer strongly 
dominated by the displacement of the Ti against the oxy- 
gen along the Ti-0 chains, but contains a significant 
component of the Pb moving against the O atoms in 
the Pb-0 planes (see Table y). Unstable Ti-dominated 
modes, similar to those in BaTi03, can be identified in 
the vicinity of the M-X line (M3/, X5 modes). How- 
ever, Pb now plays an active role in the character of the 
majority of the unstable branches, notably those termi- 
nating at M5/ and X5/ . Also, the Pb-dominated branch 
emanating from the ferroelectric T mode towards R has 
a much weaker dispersion than the corresponding, Ti- 
dominated, branch in BaTiC>3. In consequence, the un- 
stable localized ferroelectric distortion in real space is 
nearly isotropic, in contrast to the pronounced anisotropy 
in BaTiC>3. Finally, there is an antifcrrodistortive insta- 
bility at the R-point (R25 mode). As similarly observed 
in SrTiC>3 this instability is confined to quasi-one- 
dimensional "tubes" of reciprocal space running along 



the edges of the simple cubic Brillouin zone (R25 and M3 
modes and the branch connecting them). The branches 
emanating from this region stabilize rapidly away from 
the Brillouin zone edge towards, in particular, and 
X3. In real space, this instability appears as a cooper- 
ative rotation of oxygen octahedra, with strong corre- 
lations in the plane perperpendicular to the axis of ro- 
tation, and little correlation between rotations in differ- 
ent planes. The lack of interplane correlation, arising 
from the flatness of the R25-M3 branch, suggests the ab- 
sence of coupling between the oxygen motion in different 
planes. This will be discussed further in the next section. 

The ground state of PbZrC>3 is an antiferroelectric with 
8 formula units per unit cell, obtained by freezing in 
a set of coupled modes, most importantly modes at R 



and S(iiQ) |36|. The phonon dispersion correspond- 



ingly shows even more pronounced and complex insta- 
bilities than for PbTiC>3. Overall, the unstable branches 
are dominated by Pb and O displacements, with no sig- 
nificant Zr character. There is still a polar instability at 
the T point but the eigenvector (see Table |J) is clearly 
dominated by the displacement of lead against the oxy- 
gens while the Zr atom now moves with these oxygens. 
In fact, the modes where the Zr is displaced against the 
oxygens (r LO at 160 cm" 1 , M 3 ,, X 5 modes) are now all 
stable. The octahedral rotation branch is again remark- 
ably flat and is significantly more unstable at R25 and 
M3 than in PbTiOs. The antifcrrodistortive instability 
retains some one-dimensional character but spreads into 
a larger region of reciprocal space : the T25 and X3 trans- 
verse oxygen motions, related to the R25 mode, are still 
stable but with a relatively low frequency. We note fi- 
nally that the stiffest longitudinal and tranverse oxygen 
branches have been shifted to higher energy relative to 
the titanates. 



V. INTERATOMIC FORCE CONSTANTS 



In the previous section, comparisons between the three 
compounds were made by analyzing phonon disper- 
sion relations along high-symmetry lines in reciprocal 
space. A complementary, highly instructive picture of 
the quadratic-order structural energetics of the system is 
provided by direct examination of the real-space inter- 
atomic force constants (IFC). 

The interatomic force constants (IFC) are generated 
in the construction of the phonon dispersion relations; 
their computation has been described in Section II. Our 
convention is that the IFC matrix C a .p(lK, I'k') which re- 
lates the force F a (lK) on atom k in cell I and the displace- 
ment A773 (/'«/) of atom k' in cell /' is defined through the 
following expression: F a (ln) — —C a .p(lK,l'K').ATp(l'K,'). 
Moreover, the total IFC can be decomposed into a dipole- 
dipole part (DD) and a short-range part (SR) , follow- 
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ing Refs. Jl9| , ^7t . Such a decomposition is somewhat 
arbitrary but is useful for understanding the microscopic 
origin of the trends among different compounds. For con- 



venience, the atoms are labeled according to Table HI, 
as illustrated in Fig. ^. The interatomic force constants 
are reported either in cartesian coordinates or in terms 
of their longitudinal (||) and transverse (_L) contributions 
along the line connecting the two atoms. The results for 
BaTi0 3 , PbTi0 3 and PbZr0 3 are presented in Tables 
[fv|, [V]and[vi. 

First, we examine the "self- force constant," which 
specifies the force on a single isolated atom at a unit 
displacement from its crystalline position, all the other 



atoms remaining fixed. The values are given in Table IV 



The self-force constants are positive for all atoms in the 
three compounds, so that all three are stable against iso- 
lated atomic displacements. Therefore, it is only the co- 
operative motion of different atoms that can decrease the 
energy of the crystal and generate an instability, such as 
is observed in the phonon dispersion relations presented 
in the previous Section. The analysis of the IFCs will 
help us to identify the energetically favorable coupling in 
the displacements and elucidate the origin of the unstable 
phonon branches. 

Next, we discuss the ferroelectric instability at T, and 
the phonon branches which emanate from it. In barium 
titanate, it was found that the unstable eigenvector is 
dominated by Ti displacement along the Ti-O-Ti chain. 
If we consider the simple case where only Ti atoms are 
allowed to displace, we find that the destabilizing con- 
tribution from the Tio— Tii || interaction itself is nearly 
enough to compensate the Ti self-force constant (Table 
^). In addition, the fact that the Tio-Tii _L interaction is 
comparatively small can account directly for the charac- 
teristic flat dispersion along T-X and T-M and the strong 
stiffening along T-R, associated with the chain-like na- 
ture of the instability. For the true eigenvector, another 
important, though relatively small, destabilizing contri- 
bution comes from the cooperative displacement of the 
Oi atoms against the titaniums along the Ti-0 chains. 
This, together with the total contribution of the rest of 
the IFCs, is responsible for the actual instability of the 
ferroelectric Ti-dominated branches in BaTi0 3 . 

For lead titanate, the energetics of the Ti-only dis- 
placements, dominated by the Ti self-force constant and 
the Tio-Tii || and _L interactions, are remarkably simi- 
lar to those in BaTi0 3 (Table 0). However, in PbTi0 3 
there is also an important destabilization associated with 
pure Pb displacements |jq]. This can be fully attributed 
to the large difference in the Ba and Pb self-force con- 
stants, while the Ao-Ai || and _L interactions are very 
similar in the two compounds. Also, the Ao-Bo || and _L 
cation interactions are of the same order of magnitude as 
in BaTi0 3 and combine to give a surprisingly small xx 
coupling. At r, symmetry considerations permit the mix- 
ing of Ti-0 and Pb-0 displacements and in the phonon 



branches which emanate from it, thus accounting for the 
nature of the ferroelectric eigenvector. However, at X, 
M and R symmetry labels distinguish the Ti-dominated 
(X 5 , M 3 / and R25' ) and Pb-dominated (X 5 < , M 2 < and R15) 
modes, which can be readily identified in the calculated 
phonon dispersion. Also, the Pbo-Pbi coupling is much 
smaller in magnitude than the Tio-Tii coupling, which 
accounts for the relatively weak dispersion of the Pb- 
dominated branch from T to R. In the true eigenvectors, 
these instabilities are further reinforced by displacements 
of the oxygens. While the longitudinal IFC between Bao 
and Oi was very small in BaTi0 3 , there is a significant 
destabilizing interaction between Pb and Oi in PbTi0 3 , 
which further promotes the involvement of Pb in the un- 
stable phonon branches. We note that the Ti -Oi longi- 
tudinal interaction is repulsive in PbTi0 3 , but it is even 
smaller in amplitude than in BaTi0 3 and its stabilizing 
effect is compensated by the transverse coupling between 
Pb and Oi. 

In lead zirconate, the unstable eigenvector at T is 
strongly dominated by Pb-0 motion, with little involve- 
ment of Zr. This can be understood by comparing, in 
Table |y|, the energetics of Zr-only displacements with 
those of Ti-only displacements in PbTi0 3 and BaTi0 3 : 
the Zr self-force constant is significantly larger and the 
Zro-Zri || and _L interactions are smaller, so that Zr can- 
not move as easily as Ti. Also, the Zro-Oi || interaction 
is now significantly repulsive, explaining why the Zr atom 
does not move against the oxygens, but with them. As 
for the titanates, we note finally that the Zr atoms are 
mainly coupled along the B-0 chains, so that the char- 
acteristic dispersion of the B-atom modes is preserved, 
only at higher frequencies. On the other hand, the Pb 
self-force constant is much smaller, the Pbo-Pbi || and _L 
interactions are only slightly smaller, and the destabiliz- 
ing coupling between lead and oxygen is similar to that 
in PbTi0 3 , accounting for the involvement of Pb in the 
instability. 

Finally, we discuss the antiferrodistortive instability 
identified with the R25 and M 3 modes and the branch 
along R-M connecting them. There is a marked variation 
in the frequency of the R25 mode in the three compounds, 
ranging from the lack of any instability in BaTi0 3 , to 
PbTi0 3 with an unstable R25 mode that nonetheless does 
not contribute to the ground state, and finally to PbZr0 3 
in which the R25 mode is even more unstable and con- 
tributes significantly to the observed ground state | ^6| . 
The eigenvector of this mode is completely determined 
by symmetry and corresponds to a coupled rotation of 
the corner-connected oxygen octahedra. Its frequency 
depends only on the oxygen IFCs, predominantly the 
self-force constant and the off-diagonal coupling between 
nearest neighbor oxygen atoms. In fact, the latter (for 
example, Oi y -02 Z in Table VI) is remarkably similar in 
all three compounds. The trend is therefore associated 
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with the rapid decrease in the transverse O self-force con- 
stant from BaTi0 3 to PbTi0 3 to PbZr0 3 and the result- 
ing compensation of the contribution from the self-force 
constant by the destabilizing contribution from the off- 
diagonal coupling. 

The self-force constant can be written as a sum 
over interatomic force constants, according to the re- 
quirement of translational invariance: C a ^{lK,ln) — 
~Y1i>k' Ca,/3(' K ) /'«')■ It is therefore of interest to iden- 
tify which interatomic force constants are responsible for 
the trend in the transverse oxygen self- force constant. 
The suggestion that the trend is due to covalency-induced 
changes in the Pb-0 interactions can be directly inves- 
tigated through a "computer experiment." Everything 
else being equal, we artificially replace the IFC between 
Ao and Oi atoms in BaTi0 3 by its value in PbTi0 3 , con- 
sequently modifying the self-force constant on A and O 
atoms. For this hypothetical material, the A-atom dom- 
inated modes are shifted to lower frequencies while the 
frequency of the R25 mode is lowered to 40i cm -1 . If we 
introduce the stronger Ao— Oi interaction of PbZr0 3 , we 
obtain an even larger R25 instability of 103i cm -1 . 

The previous simulation demonstrates the crucial role 
played by the lead-oxygen interaction in generating the 
AFD instability. However, this change alone is not suf- 
ficient to reproduce the flatness of the R25-M 3 branch, 
as the corresponding frequencies of the M 3 mode in the 
two hypothetical cases above are 92 cm -1 and 25« cm -1 , 
respectively. Naively, the absence of dispersion of the 
antifcrrodistortive mode along that line would be inter- 
preted as the absence of coupling between the oxygens 
in the different planes. However, as can be seen in Table 
VI, the yy transverse coupling between Oi and 3 is far 
from negligible, and acts to amplify the AFD instability 
at R with respect to M. In the lead compounds, however, 
this is compensated by another yz coupling, between Oi 
and O5. The latter is significantly smaller in BaTi0 3 (by 
35 %). If we consider a third hypothetical compound in 
which this coupling in BaTi0 3 is additionally changed 
to its value in PbTi0 3 , we recover a fiat behavior along 
the R-M line. In the lead perovskites, the flatness of 
this band appears therefore as a consequence of a com- 
pensation between different interplane interactions, and 
cannot be attributed to complete independence of oxygen 
motions in the different planes. 



First, we remark that B-0 interactions depend 
strongly on the B atom, being similar in PbTi0 3 and 
BaTi0 3 , and quite different in PbZr0 3 . In fact, the SR 
force contribution to the Zro-Oi interaction and Tio-Oi 
are very similar, so that the difference arises from the 
dipolar contribution. In PbZr0 3 , this contribution is re- 
duced in consequence of the lower values of the Born 
effective charges (see Table [v]). This trend provides an- 
other example of the very delicate nature of the compen- 
sation between SR and DD forces, previously pointed out 
for BaTi0 3 [|9]j37|. 

Next, we remark that A-0 interactions depend 
strongly on the A atom, being similar in PbTi0 3 and 
PbZr0 3 , and quite different in BaTi0 3 . This change 
originates in the covalent character of the bonding be- 
tween Pb and O, which results both in smaller A-0 SR 
coupling and a larger Born effective charge for Pb. Even 
though the impact of the latter on destabilizing the DD 
interaction is partly compensated by the increased Cqq, 
the net effect is to promote the Pb-0 instability. 

As discussed above, the self-force constant can be writ- 
ten as a sum over interatomic force constants. It can be 
easily verified that the trends in the self-force constants 
observed in Table III are primarily associated with the 
trends in A-0 and B-0 interactions. 

The rest of the IFCs given in Table VI are actually 
remarkably similar. For example, A-B interactions are 
apparently insensitive to the identity of A (Ba, Pb) or 
B (Ti, Zr). This is true also for A-A, B-B and most 
0-0 interactions. The small differences observed can at 
least in part be attributed to differences in the lattice 
constants and in for the three compounds. 

The similarities in IFCs among compounds with re- 
lated compositions offer an intriguing opportunity for the 
modelling of the lattice dynamics of solid solutions. In 
the simplest case, the lattice dynamics of ordered su- 
percells of compounds such as PZT could be obtained 
by using the appropriate A-0 and B-0 couplings from 
the pure compounds and averaged values for the A-B, 
A-A, B-B and 0-0 interactions. In a more sophisti- 
cated treatment, the dipolar contributions could be sep- 
arately handled within an approach correctly treating the 
local fields and Born effective charges. Implementation 
of these ideas is in progress. 



VII. CONCLUSIONS 



VI. DISCUSSION 

In Section IV, we observed marked differences between 
the phonon dispersion relations and eigenvectors in the 
three related compounds. Through the real-space analy- 
sis in the previous section, we have seen that these differ- 
ences arise from changes in a few key interatomic force 
constants. 



In this paper, we have described in detail the first- 
principles phonon dispersion relations and real-space in- 
teratomic force constants of cubic PbTi0 3 and PbZr0 3 
and compared them with results previously obtained for 
BaTi0 3 . The modifications induced by the substitution 
of Ba by Pb and of Ti by Zr are seen to be most easily 
understood by considering the real-space IFCs. The re- 
placement of Ba by Pb strongly strengthens the A-0 cou- 
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pling, which is directly responsible for both the involve- 
ment of Pb in the ferroelectric eigenvector and the ap- 
pearance of the antiferrodistortive instability. The two- 
dimensional real-space character of the latter results from 
an additional slight modification of the 0-0 coupling. 
The replacement of Ti in PbTi03 by Zr strongly modifies 
the B-0 interaction, suppressing the involvement of Zr in 
the unstable modes of PbZr03. The decrease of the Born 
effective charges along the B-0 bonds is a crucial factor 
in modifying this interaction. In addition, the substitu- 
tion of Ti by Zr slightly strengthens the Pb-0 coupling. 
Apart from these modifications, the other IFCs are re- 
markably similar in the three compounds studied. The 
consequent prospects for transferability to solid solutions 
were discussed. 
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FIG. 1. Calculated phonon dispersion relations of BaTiOa, PbTiOs and PbZrOs along various high-symmetry lines in the 
simple cubic Brillouin zone. A color has been assigned to each point based on the contribution of each kind of atom to the 
associated dynamical matrix eigenvector (red for the A atom, green for the B atom, and blue for the oxygens) J3f|. Symmetry 
labels follow the convention of Ref. pc| , with the A atom at the origin. 

FIG. 2. Schematic three-dimensional view of the atoms labeled in Table |lll|. 



TABLE I. Lattice parameter (A) at which DFPT calculations were performed, Born effective charges (| 
dielectric constant in the cubic perovskite structure for the three ABO3 compounds. 



and optical 





BaTiOa 


PbTiOs 




PbZrOs 




4.00 


3.97 




4.12 


Z*A 


+2.74 


+3.87 




+3.93 


z* B 


+7.32 


+7.04 




+5.89 


Z o ± 


-2.14 


-2.57 




-2.50 


Z *°\\ 


-5.78 


-5.76 




-4.82 


£00 


6.75 


8.24 




6.97 


TABLE II. 


Normalized dynamical matrix eigenvector for the unstable ferroelectric mode at V (z-polarization). 


The cor- 


responding eij 


^endisplacement in real space can be obtained by dividing each value by the appropriate mass factor v ' Mi OU 


ABO3 


A B 


0, 


O, 


o z 


BaTiOs 


+0.0178 +0.6631 


-0.2842 


-0.2842 


-0.6311 


PbTiOs 


+0.2314 +0.4024 


-0.4792 


-0.4792 


-0.5704 


PbZrOs 


+0.5033 -0.1786 


-0.5738 


-0.5738 


-0.2374 
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TABLE III. Label assigned to various atoms in terms of their position in reduced coordinates. 



A 


( 0.0, 0.0, 0.0) B 


( 0.5, 0.5, 0.5) 




( 0.5, 0.5, 0.0) 


Ar 


( 0.0, 0.0, 1.0) Bi 


( 1.5, 0.5, 0.5) 


o 2 


( 0.5, 0.0, 0.5) 








o 3 


(-0.5, 0.5, 0.0) 








o 4 


( 0.5, 0.5,-1.0) 








Ok 


f-0 5 




TABLE IV. Self-force constant 


(Ha/Bohr 2 ) on the different atoms 


in the unit cell. 




Atom 


Direction 


BaTi0 3 


PbTi0 3 


PbZr0 3 


Ao 


x=y=z 


+0.0806 


+0.0247 


+0.0129 


Bq 


x=y=z 


+0.1522 


+0.1393 


+0.2302 


Oi 


x=y 


+0.0681 


+0.0451 


+0.0166 




z 


+0.1274 


+0.1518 


+0.2758 



TABLE V. Selected longitudinal (||), transverse (_L) and cartesian (a/3) interatomic force constants (Ha/Bohr 2 ) between 
different pairs of atoms. The dipole-dipole (DD) and remaining short-range (SR) contribution, have been separated following 
the scheme described in Ref. [Hi. 



BaTiOa PbTiOs PbZrOs 



Atom 




Total 


DD 


SR 


Total 


DD 


SR 


Total 


DD 


SR 


B - 


•Oi 


(II) 


+0.0094 


+0.2325 


-0.2231 


-0.0012 


+0.1865 


-0.1877 


-0.0687 


+0.1380 


-0.2067 






(-L) 


-0.0211 


-0.0430 


+0.0218 


-0.0178 


-0.0417 


+0.0239 


-0.0100 


-0.0358 


+0.0258 


B - 


-Bi 


(II) 


-0.0672 


-0.0368 


-0.0304 


-0.0615 


-0.0285 


-0.0330 


-0.0499 


-0.0211 


-0.0288 






(-L) 


+0.0075 


+0.0184 


-0.0109 


+0.0065 


+0.0142 


-0.0077 


+0.0054 


+0.0105 


-0.0052 


B - 


-0 4 


(II) 


+0.0156 


+0.0086 


+0.0070 


+0.0135 


+0.0069 


+0.0066 


+0.0106 


+0.0051 


+0.0055 






(-L) 


+0.0009 


-0.0016 


+0.0007 


+0.0015 


-0.0015 


+0.0006 


+0.0012 


-0.0013 


+0.0002 


B - 


-Ao 


(II) 


-0.0286 


-0.0212 


-0.0074 


-0.0277 


-0.0241 


-0.0036 


-0.0271 


-0.0216 


-0.0054 






(-L) 


+0.0134 


+0.0106 


+0.0028 


+0.0157 


+0.0121 


+0.0036 


+0.0145 


+0.0108 


+0.0037 






(xx) 


-0.0006 


+0.0000 


-0.0006 


+0.0012 


+0.0000 


+0.0012 


+0.0007 


+0.0000 


+0.0007 


A - 


-Oi 


(II) 


-0.0004 


+0.0114 


-0.0118 


+0.0108 


+0.0162 


-0.0054 


+0.0139 


+0.0169 


-0.0030 






(zz) 


-0.0108 


-0.0154 


+0.0045 


-0.0110 


-0.0181 


+0.0071 


-0.0103 


-0.0163 


+0.0060 


Ao- 


-Ai 


(II) 


-0.0112 


-0.0052 


-0.0060 


-0.0108 


-0.0086 


-0.0022 


-0.0094 


-0.0093 


-0.0001 






(-L) 


+0.0038 


+0.0025 


+0.0012 


+0.0054 


+0.0043 


+0.0011 


+0.0056 


+0.0047 


+0.0009 



8 



TABLE VI. Interatomic force constant matrix in cartesian coordinates (Ha/Bohr 2 ) between various pairs of oxygen atoms. 
Lines and columns of the matrix correspond respectively to x, y and z displacements for the first and second atom mentioned 
in the first column of the Table. 

Atoms BaTiOa PbTiOa PbZrOa 



O1-O2 



O1-O3 



O1-O4 



O1-O5 



+0 


.0037 





,0000 





.0000 





.0000 


— 0. 


,0087 


+0 


.0119 





,0000 


+0. 


,0274 


-0 


.0087 


-0 


.0019 





,0000 





.0000 





.0000 


+0 


,0017 





.0000 





,0000 


0. 


,0000 


+0 


.0091 


-0 


.0003 





,0000 





.0000 





,0000 


-0. 


.0003 





.0000 





.0000 





0000 


-0 


.0321 


-0 


.0006 


-0 


0013 


+0 


.0007 


-0 


.0007 


+0. 


.0013 


+0 


.0007 


+0 


.0013 


+0 


.0025 


+0 


.0013 



+0 


,0035 





.0000 


0. 


.0000 


0. 


.0000 


— 


.0091 


+0 


.0123 


0. 


0000 


+0 


.0271 


-0 


.0091 


-0. 


.0012 





.0000 


0. 


.0000 





0000 


+0 


.0022 


0. 


.0000 


0. 


.0000 





.0000 


+0. 


.0079 


+0. 


.0003 





.0000 


0. 


.0000 


0. 


.0000 


+0 


.0003 


0. 


0000 





0000 





,0000 


-0. 


.0326 


-0 


0010 


-0 


.0013 


+0. 


.0010 


-0. 


.0010 


+0 


.0012 


+0. 


.0011 


+0. 


.0013 


+0 


.0022 


+0 


.0012 



+0 


.0038 


0. 


0000 





,0000 





.0000 


— 0. 


.0065 


+0. 


,0110 





.0000 


+0 


.0229 


-0. 


,0065 


-0 


.0012 


0, 


0000 





,0000 


0. 


.0000 


+0. 


.0021 





,0000 


0. 


.0000 


0, 


0000 


+0. 


.0055 


-0 


.0010 


0, 


0000 





.0000 





.0000 


-0. 


.0010 


0. 


0000 


0. 


.0000 


0. 


0000 


-0 


0362 


-0 


.0010 


-0. 


.0013 


+0 


0010 


-0 


.0010 


+0. 


.0011 


+0. 


.0010 


+0 


.0013 


+0 


.0018 


+0 


.0011 
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